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Out-of-equilibrium evolution of quantum fields in the hybrid model with quantum 

back reaction 
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The hybrid model with a scalar "inflaton" field coupled to a "Higgs" field with a broken symmetry 
potential is one of the promising models for inflation and (p)reheating after inflation. We consider 
the nonequilibrium evolution of the quantum fields of this model with quantum back reaction in 
the Hartree approximation, in particular the transition of the Higgs field from the metastable "false 
vacuum" to the broken symmetry phase. We have performed the renormalization of the equations 
of motion, of the gap equations and of the energy density, using dimensional regularization. We 
study the influence of the back reaction on the evolution of the classical fields and of the quantum 
fluctuations. We observe that back reaction plays an important role over a wide range of parameters. 
Some implications of our investigation for the preheating stage after cosmic inflation are presented. 

PACS numbers: 98.80.Cq, ll.30.Qc 



I. INTRODUCTION 

The hybrid model has been proposed by Linde Q, Q 
as a possible inflationary scenario, several variants of the 
model have been discussed recently ,4, jSc JL; 0i ^1 ■ 
Here we consider the model in the context of preheating 
after inflation, and not inflation itself. This period has 
many interesting aspects and has received a wide atten- 

tion BiHiiiaEilliiiaiillillillillillii. 

For this period it is a reasonable approximation to ne- 
glect the coupling to gravity and the expansion of the 
universe. In hybrid inflation the inflaton transfers at the 
end of the slow roll period its energy to a Higgs field, 
whose classical potential changes from a symmetric po- 
tential to a double well potential, thereby going through 
a period of abundant particle production. In the hybrid 
model two principal mechanisms for particle production 
are at work: spinodal instability and parametric reso- 
nance. Their relative importance depends on the param- 
eters of the particular scenario. The hybrid model may 
arise naturallyin the context of supersymmetry and su- 
pergravity [EM IE ^| i it could also display essential 
features of an electroweak reheating (12L (ISj. Another 
aspect of the type of model investigated here may be 
that it simulates a first order phase transition, where the 
change in the effective potential does not arise from a 
decrease in temperature but is mediated by an effective 
field. It thereby replaces models 0, 0, [23, 0| where 
a rapid decrease of the temperature is simulated by an 
instantaneous quench. 

In this paper we address several new topics: 
(i) we include the back reaction of the produced Higgs 
and inflaton quanta onto the classical Higgs and inflaton 
fields and onto themselves in the Hartree approximation. 
This feature limits the exponential growth of fiuctuations 
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due to negative squared masses and to parametric reso- 
nance, which is unhibited if this back reaction is not in- 
cluded. In previous investigations quantum fluctuations 
have been included without back reaction i&| , with back 
reaction in the one-loop approximation [lOj and in the 
Hartree approximation with an UV cutoff in a super- 
symmetric model including axions p^ . The model in- 
vestigated there is rather special, it incorporates several 
scalar degrees of freedom of a non-minimal supersym- 
metric standard model. So on the one hand it has more 
scalar fields than the simple hybrid model, on the other 
hand the couplings are constrained. 

(ii) on the formal level we present a consistent renor- 
malization of the gap equations arising in the Hartree 
approximation, making use of the two-particle point- 
irreducible (2PPI) formalism [H HI llj. This is a 
first example of renormalization for coupled channel sys- 
tem out of equilibrium with Hartree- type self-interaction, 
and goes beyond the one-loop renormalization of Ref. 
|lO| and the unrcnormalized Hartree approximation of 
Ref. [13. 

(iii) we discuss, as previous authors, the approach to 
spontaneous symmetry breaking; as in Ref. we find 
that the classical Higgs field approaches the broken sym- 
metry minimum with an exponential behavior. However, 
we find that for sufficiently high excitations the system 
does not display a broken symmetry phase at late times; 
for such excitations the broken symmetry minima are 
washed out by quantum fluctuations. 

We are not able to discuss thermalization, here. Ther- 
malization is not supposed to occur in the Hartree ap- 
proximation when the fields are homogeneous in space. 
One of the possibilities in the present context has been 
discussed in Refs. 0, 0|: a transition between a quan- 
tum description at early times, until the the low energy 
modes reach high occupation numbers, and appending a 
classical evolution of this system at late times. For such 
a classical evolution in a <I>^ model in 1 -I- 1 dimensions 
thermalization has recently been reported and systemat- 
ically investigated in Ref. In a 3 -I- 1 dimensional 
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massless model the system is found to evolve towards 
a turbulent behavior at late times |26j . Thermalization 
and/or "equilibration" has also been discussed within a 
purely quantum field theoretical description for spatially 
inhomogeneous fields in the Hartree approximation . 27j 
and for spatially homogeneous fields, when going to the 
next-to-leading order in a ^[N expansion within the 2PI 
formahsm (2Pl-NL0-f /N) [2l[29] or to the bare vertex 
[30| approximations. The issue is not closed at present, 
and the different approaches have to be studied in con- 
junction in the future. 

The plan of the paper is as follows. We start by briefly 
reviewing the hybrid model and its application to the 
preheating stage after cosmic inflation in Sec. ^ In 
Sec. mil we present an effective action and the deriva- 
tion of the equations of motion for the nonequilibrium 
quantum dynamics in the hybrid model. In Sec. IIVI we 
present approximations to the effective action and spec- 
ify the Hartree approximation by giving the expressions 
for the fluctuation integrals and the energy contributions. 
The renormalization is discussed in Sec. In Sec. IVII 
we give details of the numerical implementation. The nu- 
merical results for different parameter sets are presented 
and discussed in Sec. IVIII We end with conclusions in 
Sec. I Villi The paper is completed by three appendixes. 



II. THE HYBRID MODEL 

A. General form of the potential 

We consider the hybrid model as proposed by Linde 
defined by the Lagrange density 
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with the potential (hybrid potential) 



A 



(2.1) 



(2.2) 



We will refer to as the inflaton and to X as the Higgs or 
the symmetry breaking field. For simplicity the depen- 
dence of all flelds on the space-time arguments x = (t, x) 
is suppressed. The discussion will be restricted to a 
space-time with Minkowski metric. The Higgs field has 
a double-well potential with a classical vacuum expecta- 
tion value given by v. Both fields $ and X are assumed 
to have non-vanishing classical expectation values (order 
parameters of the fields), i.e. 



(2.3) 
(2.4) 



The generic shape of the classical hybrid potential 
V{(j), x) is shown in Fig. Classically there are the 
two degenerate minima at the points with cj) — and 
X = ±Au. They can easily be identified in the contour 



lines at the base of the plot in Fig. ^ However, if the 
corrections from quantum fiuctuations around the classi- 
cal solutions are included, the minima of the interacting 
theory are not equal to the ones in the classical theory. 




FIG. 1: Generic shape of a hybrid potential V{(j),x) ^ ^ 
function of the classical fields (j) and x', at the base of this plot 
the contour lines of the potential are show. 

The Higgs field has an effective mass square (<I>) = 
—Xv^ + g^^^. Thus the field X, i.e. the classical part 
X and its quantum fluctuations, becomes unstable if the 
absolute value of the inflaton field <& is lower than 



(2.5) 



The region with — <I>c < $ < <I>c in the potential is called 
the spinodal region. 



B. Energy scales and parameters 

If the hybrid model is taken as an effective model for 
the (p)reheating stage at the end of cosmic inflation, the 
mass parameters and and the coupling constants 
g and A are constrained by the observations of the vari- 
ous cosmic microwave background experiments (see, e.g., 
Refs. 1^ |3^ [s^, 13 for recent analyzes). In order to 
fulfill the slow-roll condition the inflaton mass m has to 
be chosen very small Q 



(2.6) 



This means that the potential in the direction close to 
X ~ is very flat. Indeed when treating the period after 
inflation the inflaton mass can be neglected (see, e.g., 
|35l|): we have chosen = in the numerical simulations 
throughout. As discussed in Q the couplings A and g'^ 
can vary over a wide range, depending on the speciflc 
preheating scenario. This is interrelated with the choice 
of the mass scale v which is chosen to be in the numerical 
simulations. As we do not couple to gravity its absolute 
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physical value is irrelevant here. Of course this choice 
enters into the time, momentum and energy scales. 

As a concrete phenomenological model it can be em- 
bedded at least at two manifestly different energy scales. 
If one has electroweak preheating |^ 0, 0, 
mind, the symmetry breaking field mimics the standard 
Higgs sector of the standard model, while ignoring contri- 
butions from the gauge and fermion fields. Of course we 
would have to generalize X to be a complex doublet in or- 
der to represent the SU(2) symmetry breaking Higgs sec- 
tor. The vacuum expectation value v in this case would 
be chosen equal to 246GeV. The classical Higgs mass is 
= Xv"^. If the phase transition at the end of infla- 
tion takes place at the scale of a Grand Unified Theory 
(GUT), the symmetry breaking field X could be, e.g., a 
very heavy sneutrino, i.e. the scalar super-partner of one 
of the light neutrinos. The reheating in such a scenario 
would come from the decay of this heavy sneutrino. We 
will not restrict our study to one of these two scenarios 
in the following. 



III. EFFECTIVE ACTION 

A. Quantum fluctuations and their 
non-perturbative resummation 

The time evolution of out-of-equilibrium systems in 
quantum field theory is described within the closed-time- 
path (CTP) or Schwinger-Keldysh formalism [s^ Is^ . 
Nonequilibrium quantum field theory has seen a tremen- 
dous progress, after seminal pubHcations in the 1980's 
m, mainly in the last decade [IS IH El El El El ■ 
The main ingredient, besides the CTP formalism itself, 
are various approximations based on variational princi- 
ples. The one-loop, Hartree and and large- TV approxima- 
tions can be based on the two-particle irreducible (2PI) 
formalism E3jE3- This formalism has also been used be- 
yond these mean field approaches ,283 30j J£, . As 
an alternative approach beyond leading orders we have 
recently adapted '23', '23| the so-called two-particle point 
irreducible (2PPI) formalism 22, 52] to nonequilibrium 
quantum field theory. 

Within the framework outlined above there are sev- 
eral ways for taking into account the quantum fluctua- 
tions and their back reaction onto the classical (mean) 
fields and onto themselves. In any case they should to 
be treated non-perturbatively, as the amplitudes of the 
fields can become very large. In practice one resums an 
infinite number of perturbative Feynman diagrams via a 
Schwinger-Dyson equation. In this section we will derive 
all equations of motion, i.e. the classical and quantum 
equations, from an effective action that incorporates full 
back reaction of the quantum modes. 

Within nonequilibrium quantum field theory the ex- 
istence of a conserved energy is automatically guaran- 
teed if all equations of motion are derived from the same 
variational functional, i.e. the effective action. In the 



context of relativistic quantum field theory Cornwall, 
Jackiw and Tomboulis (CJT) have derived the 2PI ef- 
fective action in Ref. |43|. In addition to a local source 
term J{x) for a mean field 4'{x) one introduces a bilo- 
cal source term K{x^x') for the Operator {^{x)^{x')) — 
G{x, x') + (t>{x)(j){x') in the generating functional. The ef- 
fective action F follows from a Legendre transformation 
of the generating functional with variational parameters 
0(2;) and G{x,x'). The leading orders in the 2PI or CJT 
action exactly reproduce the well know one-loop, large- iV 
and Hartree approximations. 

A related effective action is the 2PPI effective action 
IHIs^. There one has a local source term K{x,x') = 
K{x)5{x — x') for the composite operator ($^(a;)) thus 
making all self-energy insertions into the propagator 
G{x, x') local. In the 2PPI effective action formalism 
the renormalization of the Hartree (and leading order 
large- iV) approximation is very transparent, because the 
non-perturbative counterterms in the effective theory are 
mapped one-to-one to the standard counterterms of a 
perturbative expansion of the basic Lagrangian. In addi- 
tion the 2PPI scheme allows a systematic improvement of 
the Hartree approximation without all the computational 
complications from nonlocal self-energies as introduced 
beyond leading oder in the 2PI scheme. 



B. The two-particle point-irreducible formalism 

At first, a slightly more general Lagrange density for a 
N component field $ — ($1, . . . , ^n) is given by 



t^i^) = \d^^\x)d^''i>'(x)-V{x) , 



(3.7) 



V{x) = im?-$*(x)$^(a;) 
2 ■' 



+ j^ X,,ki^\x)^^ {x)^\x)^\x) , (3.8) 

where a summation over i, j, k,l = I, . . . , N etc. is under- 
stood and we have introduced a coupling constant ma- 
trix Xijki and a mass matrix mfj. The fields may have 
classical expectation values, i.e. = 

The standard hybrid potential in Eq. H2.2() follows from 
Eq. H3.8|l for the case N = 2 with the identifications 

(3.9) 
(3.10) 
(3.11) 
(3.12) 
(3.13) 
(3.14) 
(3.15) 

= A2221 = . (3.16) 

The 2PPI formalism is formulated in terms of the mean 
fields (j>^ and local propagator insertions Aij. The lo- 
cal insertions in the propagator are resummed via a 
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Schwinger-Dyson or gap equation. The Green's function 
G fulfills the local equation 

{G-%{x,x')=i{n5^j+Ml{x))5^''\x~x') , (3.17) 

where M-fj is a variational mass parameter explained be- 
low. 

For the general potential in Eq. H3.8|) the 2PPI effective 
action can be written as (see e.g. Ref. (s^. Is^ l 



+ i\^3kl I d''xX\x)A'''ix) (3.18) 



kl, 



with 



A'\x) = -2 



5T 



2PPI 



SMUx) ■ 



(3.19) 



The masses Mfj{x) have to fulfill the so called gap equa- 
tions 

Mlix) = ml + {!^Hx)<i}(x) + A'='(x))(3.20) 

The term r^^^^[0%A^y denotes the infinite sum over 
all 2PPI diagrams. A diagram that does not fall apart 
if two lines meeting at the same point are cut is called 
two-particle point-irreducible. 

The classical action S[(f>^] is given by 



d'^x 



-d,fix)d^^Hx) - -ml4,\xW(x) 



-Kjki4>\x)4>^{x)c^''{xW{x) 



(3.21) 



The term proportional to an integral over /V^{x)A^^{x) 
in Eq. H3.18|l corrects the double counting of bubble 
graphs. 

With the help of the identifications in Eq. (pTl^ - (prTH|) 
the 2PPI effective action for the hybrid potential in 
Eq. 12. 2|l reads explicitly 

r[(/>,x, A^^, A^;^] 



+ ^ / d^x(A^^(x)A^^(x)+2A2jx)) .(3.22) 



The gap equations for the masses M.^ are 

Ml^{x) = + {x\x) ^ A^^{x)) , (3.23) 



Ml^{x) = -\v^+2,\{x\x) + A^^{x)) 



•2 '-2,\{X^' 
+g^ {4>\x) + A^^(x)) , (3.24) 
Ml^{x) = Ml^{x) = 2g^ Wx)x{x) + A^^(a;))(3.25) 



The self-energy insertions A are 

Jp2PPI 

A^^{x) = 
A^^{x) = 



^SMl^x) 
'sMl^ix) 

Jp2PPI 



Inverting the gap equations gives 



(3.26) 
(3.27) 
(3.28) 



Aci,^{x) = 
A^^{x) = 



^[Ml,ix) + Xv--'^{MUa 

-<l>^ix) , 

l{Ml^{x)-m^)~x'{x) , 

1 
2^ 



9 



(3.29) 
(3.30) 

(3.31) 



The effective action can therefore be expressed in a form 
without the A's. 

T[<j>,x,Ml^,Ml^,Ml^] 
1 



-Im14x)^\x) g'cj,\x)x'{x) 

-lMl^i^)x'{x) + ^x\x)-^v' 
-Ml^{x)(l>{x)x{x) 



+r2PPi [<p,x.Ml^,Ml^,M 



' 1 
xxJ 



(3.32) 



The latter form of the effective action makes it easy to 
derive a renormalized effective action. 



C. Equations of motion 

The equations of motion for the classical fields </> and 
X follow from 



^0 and ^ = 0. 



S(j){x 



6xix) 



(3.33) 
(3.34) 



If we take only the causal part of all functions by restrict- 
ing them to the positive time branch of a closed time path 
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(CTP; for the Schwinger-Keldysh closed time path for- 
mahsm see Ref. [s^ y^, |3 Ea), we get the following 
classical equations of motion 

^p2PPI 







(3.35) 



ax + M'x + MLcl>--2Xx 



ST 



2PPI 



Sx 



(3.36) 



The quantum equations of motion, i.e. the Schwingcr- 
Dyson equations or gap equations, are derived from 

ST ^ ST . . ST 



0, 



SMl^ix) 



and 



0(3.37) 



They have already been denoted in Eq. H3.23|l - H3.25|l . 
In general the gap equations form a coupled system 
of self-consistent nonlinear equations. Due to the self- 
consistency an infinite number of one-particle irreducible 
(IPI) diagrams is resummed. 

In the following we will assume spatially homogeneous 
fields. The classical fields obey (p — (j){t) and x = x{t) ^^'^ 
the Green's functions can be expressed via their Fourier 
components 



Gy(t,i';x,x')= J 



jD-1 



p 



(2^) 



D-1 



)Gy(t,t';p)(3.38) 



Since the 2PPI formalism resums local self energy inser- 
tions in the Green's function, Gij(t,t']p) can be rewrit- 
ten in terms of mode functions fi{t;p) and fj{t';p). If 
we do so, we have a Wronskian matrix, that is diagonal. 
The quantization of the theory at t = is explained in 
more detail in Sec. IV Al The Green's function reads 

2 



a=l 



fr{t-p)fr{t';pmt-t') 



+ft{t';p)frit-^pMt'-t) 



(3.39) 



where the mode functions satisfy 

ft{t;p) + [p'S,,+Ml{t)] ff{t-p) = . (3.40) 

The fundamental solutions of this system of coupled dif- 
ferential equations will be labeled with Greek letters 
a, /?, . . .. So there are four different complex mode func- 
tions. 

In Eq. H3.39|l we have introduced the quantities Wq, 
defined by 



- . /r^2 



(3.41) 



They will be explained below, as they depend on the 
initial conditions via the initial masses TOq „ . 

The explicit form of the equations of motion depends 
on a given approximation for r^^^^. In addition one has 
to renormalize the equations of motion for all approxi- 
mation schemes beyond the classical, i.e. the tree level. 



IV. APPROXIMATIONS 

The infinite sum of two-particle point-irreducible dia- 
grams in r^^^^ can be truncated in several ways. We will 
discuss a loop expansion here. The loop expansion can 
be formulated as 



^2PPI 



£=1 



(4.1) 



where the index ^ = 1,2,... corresponds to the number of 
loops in a 2PPI diagram. The dots in the last equation in- 
dicates all contributions beyond the two-loop order. The 
fluctuation integrals A^- are expanded analogous as 



A.,(t) = A«(0 + Ag)(t) + 



(4.2) 



The one-loop order in the 2PPI loop expansion is equiv- 
alent to the Hartree approximation, as already stated 
above. 

a. Zero-loop - classical approximation. To zero-loop 
order one discards the term r^^^' completely, which leads 



to Aij = and ^ 



0. 



The zero-loop contribution to the energy denotes 
1 



-Mi4m'it)~g'r{t)x'{t) 



1 



id it) 



MUt)x\t) 



X 



A 



+Mi^imt)xit) 

3A 



+- 



{Ml,{t)-m^y 



(4.3) 



1 



{Ml^{t)-m^) {Ml^{t) + W) 



b. One-loop - Hartree approximation. The one-loop 
approximation is identical to what we will call Hartree 
approximation. The sum of all two-particle point- 
irreducible diagrams is truncated at r^^^^ P^^^ with 
(see the diagram in Fig. ^ 



(4.4) 



T^''>[Ml^,Ml^,Ml^] = ^TT\n [g-'] 
We have introduced 

g = C;^^ ^^x I (4 5) 



Gy^tf, G 



XX 



The functional derivative of P'^^ with respect to A^|- 



gives the tadpole insertions 
1 f d'^-^p ' 



2 

E 

a=l 



(27r)^-i 



G,j{t,t-p) + Gj,{t,t-p) 
1 



(27r 



2uJr 



(4.6) 

Re[/f(t;p)/;*(t;p)] . 

(4.7) 



6 



Q 
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FIG. 2: One-loop bubble diagram representing the Hartree 
approximation. The line denotes a propagator Q as given by 
Eq. (HSJ 



(1) 



Ag and A« = 



There a three different insertions, A 

A'^J. The gap equations H3.23|l - H3.25|l resum one-loop 
bubble graphs, because the Green's functions Gy depend 
on the variational mass parameters A^f,-. At one-loop 



order we have 



0. 



(50 Sx 

The contribution of the bubble graphs to the energy is 
defined by the relation 



6T^^nMl.,Ml^,MUdMUt) 



dt 



(27r) 



-^G.,j{t,t- p) 



dt 

dMlit) 
dt 



-(4.8) 



This equation can be integrated explicitly if one uses 
the equations of motion for the mode functions f"{t;p), 
yielding 



jD-1 



xlRe f${t;p)f^*it;p) 



+Re[f^it;p)f^*it;p) 

+ (P' + Ml4t)) Re [f^it;p)f^*it;p)] 
+ {p'+Ml^{t)) Re [f^{t;p)f^*{t;p)] 

+2Ml^me[f^it;p)f^*it;p)] \ .(4.9) 



The momentum integrations in the quantities A^^^ and 
are divergent and thus have to be renormalized prop- 
erly. We will discuss this issue in the next section. 

c. Two-loop - sunset diagrams The only two-loop 
diagram appearing is the sunset diagram displayed in 
Fig. 13 This sunset graph leads to time integrations over 
the past of classical and quantum fields ("memory in- 
tegrations]^ and introduces scattering of the quanta (see 
Ref. for two- loop simulations of theory in 1-1-1 

dimensions) . While it would be interesting to study how 
this next-to-leading order diagram affects the dynamics 
studied here, this is beyond the scope of this work. 



FIG. 3: Two-loop sunset diagram; the solid lines denote the 
propagator Q, the dashed lines the classical fields ip = {(j>,x)- 



V. RENORMALIZATION 

A. Initial conditions 

The choice of initial conditions for the quantum sys- 
tem has to be discussed together with its renormalization. 
We will take a Gaussian initial density matrix with non- 
vanishing initial values (/)(t = 0) and x(i = 0) for the 
classical field amplitudes. For the renormalization of the 
equations of motion we need a properly quantized system 
at the initial time. In order to satisfy the usual canonical 
commutation relations for the creation and annihilation 
operators of the quantum fields, we choose a Fock space 
basis at i = 0. The basic quanta are defined by diagonal- 
izing the mass matrix at t = and by choosing canonical 
initial conditions (see below) for the mode functions (see 

e.g. iiliaiiiiii) 

We define the initial masses toq.q as the eigenvalues of 
the initial mass matrix A^y (0), i.e., by the equation 



mlJt{0;p)~Ml{0)f^{Q;p) = 
The eigenvalues are given by 



(5.1) 



MlJO) + MlJO) 



± 



(5.2) 



We denote the corresponding eigenvectors by Oia , where 
the index a refers to the eigenvalue, and the latin indices 
to the components. The canonical initial conditions at 
t = for the mode functions are 

/f(0;p) = (5.3) 
/f(0;p) = -ic^„O,„ = -zc^„/f(0;p) . (5.4) 

The Wronskian matrix of these mode functions is then 
given by 



m/r,/f) = [/r*(0;p)/f(0;p) 
-/r*(0;p)/f(0;p)] 



(5.5) 
(5.6) 



As the eigenvectors are orthogonal this matrix is diago- 
nal. Choosing the normalization 



OiaOifj = 5a(} 



(5.7) 
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the Wronskian matrix becomes 



1,2 



= -2i 



ui 

LU2 



If M'^^{0) > 7W^^(0) we can fix the matrix O as 



O = 
t'dwd = 



cos I? sin -d 
— sini? C0S1? 

1 



(5.8) 



(5.9) 



■^lx(o)--^L(o) 



0) 



For the opposite case A^^^(O) < M^^{0) one should in- 
terchange ttiq I with uIq 2 and switch d — > — t?. 

B. Isolation of the divergent contributions 

The propagator insertions A,-j^ in Eq. H4.6|l are diver- 
gent. First we have to isolate the divergent contributions 
from the finite parts. Then the renormalization can be 
performed within a suitable regularization scheme. 

A strategy for the isolation of the divergences via the 
perturbative expansion of the mode functions in terms of 
partial integrations has been described in Refs. 
We have summarized what we need in Appendix IbI 

In dimensional regularization with the abbreviation 



Le = - — 7 -I- ln47r we have, using Eq. ljB8|) from Ap- 
pendixtel and Eq. HA1|I and HA4p from Appendix 1^ 



1 f d^'-'p 



2 J (27r)^-i 



G„(i,i;p)+Gj.(i,i;p) 
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1 I OiQ,OjcK 
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167r2 
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16^ 
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In- 



X OiaOjfjOipOkaVklit) 
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IGtt 



2 [L, + l]Mlit) 



^0„Oja In ■ 



'0, a 



IGtt 



X In- 



a, 13 



m, 



■In- 



/i2 ml^-mlj; m2 



0, a 



(5.11) 



where '^ia'^ja'^o, q = -^ij(O) has been used. The 
potential (t) is defined as 



H,(t) = MUt)-Ml,{0) . 



(5.12) 



We can define the finite part of Alj-* by a subtraction as 
[see Eq. (jBSj] 



,(1) 



„(i) 



(27r) 



D-1 



a,f3 



1 

2u}o 



Re[fnt;p)f--*{t;p)] 



(5.13) 

OicOjj3 Oifi Oka ^kl (t) 



Oio-Oja^aP 

1 

+ 7 T 

UJi3[UJa + LUp) 

The momentum integrations in this expression are con- 
vergent, because the subtracted terms cancel exactly the 
divergent parts. There are different finite contributions 
from the divergent part of A-j-* (see. Eq. (|5.11|l ). We 
find it useful to define the following quantities 



C — 



''0, a 



In- 



'0, a 



(5.14) 



167r2 ^ 

a,/3 



OiaOjpOkaOi 



In- 



'0, a 



0,f3 



-In- 



0,/3 



0,a 



— m\ 



0,0 



TO, 



0, a 



5.15) 



and 



pi 



(1) 



.(5.16) 



Thus the full 2PPI insertion A}.' takes a very simple 
form given by 



A«(t) 



1 



167r2 



(0)) 

{Ml^{t) - Ml^lfi)) . (5.17) 



As one can see the divergent part of A^j' is directly pro- 
portional to A^lj, i.e. with a uniform factor for all com- 
binations of i and j. In particular it is independent of 
the masses TOq q, and the matrix Oij and thereby of the 
initial conditions. On the other hand the different finite 
parts depend on the initial conditions via the constants 
and involve all components of the effective mass ma- 
trix 



C. Suitable effective counterterms for the gap 
equations 

The divergent part in Eq. H5.17|l will be removed by 
effective mass counter terms. We explain in Appendix O 
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how the effective mass counterterms we are using here are 
related to standard mass and couphng constant countert- 
erms in a standard counterterm Lagrangian. More pre- 
cisely the 2PPI formalism establishes a one-to-one map- 
ping between both counterterm approaches. 

In order to have a SM^ counterterm in the gap equa- 
tions, we add to the effective action in Eq. (|3.32|) a coun- 
terterm of the general form 



+2S^^^{Ml^f (5.18) 



With the introduced effective mass counterterms the 
renormalized gap equations take the form 



^igHC^^Ml^^it) (5.19) 

-^R,xxW = -A^' + .g^'^'W + A^^i)) 

+3\{x'{t) + A^^it)) (5.20) 
-^gH^^Ml^^t) - 12A<5exx-M|,,xxW 

-SgHC^^Ml^^it) . (5.21) 



finally given by 

Ml^4t) = m,^+g^x^{t) + A^^Mt)) 

+9^C°^ (5.23) 

VC^x (■A^R.xxW~-^R,xx(0)) 
VC^x(-^R,0xW--^R,0x(O)) ' 
■^IxxW = -Xv^+g^q^^{t) + A^^Mt)) 

+9^C% (5.24) 

+g^Cl^ (A^|,xxW--^R,xx(0)) 

+9^Cl^ (-^R,0xW--^R,0x(O)) 
+3A (x'(i) + A^^,fi„(t))+3AC°^ 

+3ACi^ {Ml^4t)-Ml^40)) 

+3AC^, {Ml^^{t)~^Ml^^iO)) 

■A^R^xW = (0(<)x(t) + A0x,fin(i)) 

+2g^C°^ (5.25) 
+2g^Cl^ {Ml^4t)-Ml^40)) 
+2g'Cl^ {Ml^^{t)-Ml^^iO)) 
+2g^Cl^{Ml^^{t) - Ml^^iO)) . 

This system of linear equations is similar to the one ap- 
pearing in the 0(iV)-model in the Hartree approxima- 
tion l5fT|. It has to be solved at each time. However, 
the coefficient matrix can be diagonalized with a time- 
independent rotation matrix, because it is time indepen- 
dent itself. Such a rotation matrix is analogous to the 
factor C = (1 + ttttIh^)"^ in the renormalization of 
the 0(A^)-model in the large- approximation 46]. 



By inserting Ay (t) from Eq. H5.17|l one can see that the 
gap equations become finite if we choose 



^^00 = S^xx = S^x = - ^ + 1] (5.22) 
= -^(7-^+1 + 1-4.) . 



This choice of corresponds to a MS prescription. In 
particular, the renormalization scheme is mass indepen- 
dent. Thus the system of renormalized gap equations is 



D. Renormalized energy 

Within the Hartree approximation the contributions to 
the energy introduce logarithmic, quadratic and quartic 
divergences. The divergences have to be compensated by 
the already fixed counterterms of the last section. 

If the effective masses are identified by the renormal- 
ized effective masses of the previous section, then the 
zero-loop contribution to the energy [see Eq. (|4.3|) ] is au- 
tomatically renormalized. 

In the following we use the rotated potential V defined 
by 

Va/3(i) = OkaVklit)Ol0 (5.26) 
= OlaOl/3V00(t) + 02a02/3Vxx(t) 

+ (Oi„02/3 + 02aOi/3)V^xW (5.27) 

According to the expansion of the mode functions in Ap- 
pendix ^ the divergent part of the one-loop contribution 
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from the bubble graphs to the quantum energy is given 
by [see Eq. (jHT^ ] 



E 



(l),div(^) - 



1 



2 J (27r)^ 



26^2 +Vaa(i) 



E 

/9 



V„/3(t)Va/3(i) 



(5.28) 



The first term is quartic divergent. Its renormalization 
corresponds to a renormalization of the cosmological con- 
stant A; it is therefore somewhat arbitrary and can be 
omitted. 

If the divergent parts are evaluated in dimensional 
regularization [using Eq. (El), (EH) and (UHl) in Ap- 
pendix ^ the full one- loop contribution in Eq. 14.9|) de- 
notes 



'0, a 



647r2 



- In ■ 



L, - In ■ 



"O, a 
1^ 



a, P 



- In ■ 



'0, a 



+ - 



m, 



0,13 



^0, /5 



In- 



^0, a 



(5.29) 



where the finite part has been defined as 

^(1).. _ 1 / d^-^p ^ 1 
fi" ^ ' 2 J (2^)^-1 ^ 2l^„ 



xjRe 

+Re 



f^it;p)fr(t'P) 



f^it;p)f^*it;p) 

{p'+Ml^it))Re[f^{t;p)f^*{t;p)] 
{p' + Ml^{t))Re[f^{t;p)f^*it;p)] 
2Ml^me[f^{t;p)f^*{t;p)] 

2w^- Vaa(i) 

1 



+E 



Vc«/3(i)Va,3(i) M5-30) 



From the definition of TOq ^ and Vq/j (i) one can prove the 
identity 



E 



2"^0,aVaa(t) - X! ^"Z? (t) 
(3 



-A^2^.(t)A^5(t) 



(5.31) 



so that the divergent part becomes very simple. The full 
quantum energy E^^') is then given by 



E'^'Ht) = El}^it)--^[L, + l]Mlit)Ml{t) 



•sr^ '"0, a 



In 



1 
2 



Vaa [t] In ■ 



'0, a 



^ 327r2 

OL 



a,/3 



In- 



0,/3 



In- 



/i^ TO, 



0,a '"-0, /3 



0,/3 



'0, Q 



(5.32) 



The divergent part of i?'-^-' is proportional to Ad'^jAd'^j — 
Ml^ + M^^ + 2M'^^, i.e. the counterterm in Eq. 
as it has been expected. 

Within the given approximation we can write down the 
renormalized total energy as 



Eu 



E^°\t) + E'^^\t) - 6M^ 



(5.33) 



647r2 



'0, a 



In- 



Q,/3 



In- 



0,/3 



'0, a 



'■0,f3 



In- 



(5.34) 



Because we have added the counterterm 6A4'^ to the ef- 
fective action F, i.e., Fr = F + SM^, it has a negative 
sign in the renormalized energy. 



E. Renormalized equations of motion 

In summary one has to solve in the Hartree approx- 
imation the following renormalized equations of motion 
numerically. The classical equations of motion are given 

by 

-2gWmt) , (5.35) 

= m+Mi^^{t)x{t)+Mi^^m{t) 

-2Ax3(t)-2<7V(0xW , (5.36) 
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while the equations for the mode functions denote exphc- 
itly 

- f^{t;p)+py^{t;p) 

+ Ml^4t)f^it;p) + Ml^^it)f^{t;p) (5.37) 
= f;{t;p)+p'f^{t;p) 

+ Ml^^{t)f;{t;p)+Ml^^{t)f^{t;p)(,5.38) 

In addition one has to solve, at each time, the 3x3 
system of renormalized gap equations (|5.23|) - (|5.25|l for 
the masses TW^,-. 



VI. NUMERICAL IMPLEMENTATION 

The masses TOq ^ (a — 1,2) and the mixing angle 
have to be determined self-consistently at the initial time 
t = 0. The renormalization scale /i for the finite parts is 
fixed to /^^ = \v^. 

The equations of motion (|5.35|I - H5.38|I are solved us- 
ing a standard fourth order Runge-Kutta algorithm. We 
use a time discretization At = 0.0003. The momen- 
tum integrations are carried out on a momentum grid 
with a non-equidistant momentum discretization. We 
use a pragmatic momentum cutoff Pmax = 12 and Hp — 
300 momenta for the convergent fluctuation integrals in 
Eq. (|5.13() and (|5.30() . The accuracy of the numerical 
computations is monitored by verifying the constancy of 
the Wronskians and of the total energy. 

VII. RESULTS 

A. General outline 

The basic conception in the Hybrid model is an efficient 
energy transfer from the inflaton to the Higgs degree of 
freedom mediated by a phase transition and the associ- 
ated spinodal regime. The Lagrangian is constructed in 
such a way that this type of behavior can be expected. 
It is then a question how these expectations are real- 
ized; as mentioned in the Introduction this has been in- 
vestigated in various approximations on the classical or 
quantum level. Here we present numerical simulations in 
the Hartree approximation, which encompasses spatially 
homogeneous classical fields, and quantum fluctuations 
with a specific back reaction among themselves. Specifi- 
cally we are interested 

1. to see on which time scale and in which form the en- 
ergy transfer between the inflaton and Higgs fields 
takes place 

2. to conclude on the structure of the effective Higgs 
potential after this energy transfer. Though the 
system does not go right away into a thermal equi- 
librium phase the behavior at intermediate and late 
times can be thought as reflecting the shape of an 



effective Higgs potential, with a symmetric or bro- 
ken symmetry structure. 

3. in the spectra of the different quantum modes re- 
flecting the mechanism of particle production 

4. in finding out to which extent the transition to a 
classical description may be justified in a certain 
momentum range. 

The answer to these questions obviously depends on 
the parameters chosen for the simulations. In order to 
study the influence of the coupling strength and the 
self-coupling A on both the classical and quantum com- 
ponents of the Higgs and the inflaton fields we have per- 
formed simulations with = 0, = 1, A = 1 fixed, 
while is equal to 2A, O.IA and O.OIA (Fig. HE] and Ej) 
and for g^ = 2A with a smaller coupling A = 0.1 (Fig.|SJ. 
We have chosen x(0) = 10~^, i.e., a very small value for 
the initial amplitude of the classical Higgs field, in order 
to trigger the "spontaneous" symmetry breaking. The 
initial amplitude for the inflaton field has been fixed for 
all cases to (/)(0) = 1.697(/>c. The energy contributions for 
the simulation in Fig. 0] are displayed in Fig. 

B. Time regimes and exponential growth 

We have identified three time regimes in the simula- 
tions that we want to investigate further in the following. 
The regimes are described as follows: 

(I) Initial period, end of slow roll. Here M^^{t) > 0. 
A phase of slow rolling of the inflaton field after the 
main period of inflation. The quantum fluctuations 
still are almost negligible. 

(II) Early times, spinodal regime, M^^{t) < or os- 
cillating several times around zero. Spinodal am- 
plification of Higgs quantum fluctuations and and 
exponential growth of x{t)- 

(HI) Intermediate and late times, Ai^^{t) > and oscil- 
lations of the classical fields. Excitation of inflaton 
and mixed quantum fluctuations, parametric reso- 
nance bands in all momentum spectra. 

The first period (I) is easy to identify in Figs.EHHl only 
the inflaton decreases with time in smooth way while the 
Higgs mean field is still practically zero. 

In the early time period (II) the inflaton field passes 
through zero once or several times, depending on the cou- 
pling g^, see Figs. 0H7| The period is identified by an 
increase of |x(^)l a-nd ends once x{t) begins to oscillate 
in a regular way. A closer analysis shows that the am- 
plitude of the classical Higgs field growths exponentially. 
In Fig. Owe display on a logarithmic scale the absolute 
value \x{t)\ for simulations with g^ = 2A, A = 1, = 0, 

= 1 and x(0) — 10^^ fixed, while the initial amplitude 
of the inflaton field, 0(0), is varied from 1.2 to 1.8. The 
exponential growth sets in when (j){t) becomes smaller 
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FIG. 5: (Color online) Energy contributions for the simulation 
with parameters in Fig. 2] the red solid line denotes the total 
energy, the green dashed line the classical energy and the blue 
dotted line the quantum energy 
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FIG. 6: (Color online) The same as Fig. but for the cou- 
pling case = 0.1 A and the initial value 0(0) = 5.366 



FIG. 4: (Color online) Time evolution for the simulation with 
= 2A. Initial values: (j>{0) = 1.2 and x(0) = 1.0 x lO"'^; 
other parameters: = 0, A = 1, i;'^ = 1; we plot as a func- 
tion of time (a) the classical fields (j}{t) / (j>c (red solid line) and 
x{t)/v (green dashed line), (b) the effective masses Mij{t) 
with ij — cj>(j> (red solid line) , ij — XX (green dashed line) and 
ij = cj>x (blue dotted line) (c) the fluctuation integrals Aij(t) 
with ij — (j>4> (red solid line) , ij — xX (green dashed line) and 
y = 0X (blue dotted line); the vertical dotted lines indicate 
the times where t is equal to 19.5, 51.3 and 100.2 

than the critical value (f>c (see Eq. H2.5|l ) and stops when 
x{t) reaches the turning point which is at ~ 1. 

There does not seem to be a systematic trend for the 
dependence of the period of growth on 0(0). 

The regime (II) can be very short. For the simulation 
with a small coupling = O.OIA the transition to the 
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FIG. 7: (Color online) The same as Fig. but for the cou- 
pling case g^ = O.OIA and the initial value 0(0) — 16.97 
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FIG. 8: (Color online) The same as Fig.|I|i but for A = 0.1 



broken symmetry phase can take place within a single 
oscillation of the inflaton field (see Fig.|7^). 

The intermediate and late time period (III) is charac- 
terized by oscillations of both the inflaton and the Higgs 
field, with essentially constant period and amplitude (see, 
e.g., Fig. ISj). The Higgs mean field may oscillate around 
a nonzero value, related to a broken symmetry minimum 
of an effective potential (see below) or around x = 0, to 
be identified with symmetry restoration. This period is 
further analyzed in the next subsection. 




FIG. 9: Time evolution of the absolute value |x(i)| for simula- 
tions with the parameters = 2A, — 0, v'^ = 1 and A = 1 
and four different initial values <j!>(0) equal to 1.2 (solid line), 
1.4 (dashed line), 1.6 (dotted line) and 1.8 (dashed-dotted 
line); the corresponding arrows pointing at the t-axis indicate 
the time when (j>{t) drops below (j>c 



C. Late time averages — phase transition 

The amplitudes of the classical fields 4> and x de- 
crease very slowly, if at all, at late times, i.e., once they 



have started to oscillate in a kind of effective potential. 
Though we make no attempt to reconstruct such a po- 
tential in detail, the oscillations allow to conclude on the 
minimum and the range of such an effective potential for 
both the Higgs and inflaton fields. In this sense we can 
speak of a symmetric or broken symmetry phase for the 
Higgs field, if the minimum of its effective potential is 
at X = and X 9^ 0, respectively, and we associate this 
minimum with the time average of the Higgs field at late 
times. The shape of the effective potential depends here 
on the energy density (in place of the temperature) and 
therefore on the initial value of the infiaton field. The 
question of spontaneous symmetry breaking and of the 
point of the phase transition reduces therefore to finding 
the the dependence of x{t 00) on the initial value 0(0). 

In order to study this issue we have performed a series 
of simulations where we have varied the initial ampli- 
tude (/)(0) while keeping all the other parameters fixed. 
In Fig.^lthe time evolution of x(^) for simulations with 
0(0) — 1.9, 2.0 and 2.1 is displayed. The other parame- 
ters are 5^ = 2A, A = l, 11 = 1, m = and x(0) = 10~^. 




100 200 300 400 500 600 700 800 
t 

FIG. 10: (Color online) Time evolution of x(i) with param- 
eters as in Fig. El but for the initial values 0(0) equal to 1.9 
(dashed line), 2.0 (green solid line) and 2.1 (dotted line). 

A first inspection suggests that there is a phase tran- 
sition between 0(0) = 2.0 and 0(0) — 2.1, as for the 
latter simulation the field x oscillates around zero. From 
the simulation with 0(0) — 2.0 in Fig. EH it becomes 
apparent that the field x can jump several times from 
one "minimum" to the other if 0(0) is close to the crit- 
ical point of the phase transition. This is typical for a 
first-order phase transition and has been observed in the 
scalar 0{N) model in the Hartree approximation as well 

In Fig. ^] we display the time averages |x(c«)l and 
M'^^xi^) a.s a function of the initial amplitude 0(0). 
The other parameters are fixed to — 2 A, X — 1, v — 1, 
m = and x(0) = 10~^. A first-order phase transition is 
signaled by a non-continuous drop of the minimum value 
|x(oo)| and the effective mass from finite (positive) 

values to zero. 
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FIG. 11: The late-time absolute value |x(oo)| (diamonds) and 
the late time effective mass M'^^{qo) (triangles) averaged at 
times t ~ 1000 as a function of the initial amplitude <^(0); 
other parameters as in Fig. 2] 



D. Momentum spectra 

Using the amplitudes /"(i,p) one may define various 
"power spectra" . One of those is the integrand of the 
fluctuation integrals Aij{t), i.e., the tadpole contribu- 
tions. We have already introduced the kernel 



G,,{t,t,p)^{^*{t,pMt,p))/V 



(7.39) 



in terms of which we define the power spectrum of the 
fluctuation amplitudes 



P 

27r2 



(7.40) 



by including the momentum phase space factor. 

In Fig. 1121 we display this spectrum for the simulation 
in Fig. El at the times t equal to 19.5, 51.3 and 100.2. 
These time steps are indicated in Fig.0]by vertical dashed 
lines. Actually we have subtracted the free field part 
and the first order perturbative part of this kernel, in 
analogy to the right hand side of Eq. 15.1311 . The free 
field part rises linearly with momentum; for p — 2, the 
maximal value used in our plots, it is has typical values of 
0.05 and would be visible. It makes the tadpole integrals 
divergent; as discussed above, in our computations this 
divergence is absorbed by dimensional regularization and 
renormalization. 

At early and intermediate times the Higgs fluctuations 
dominate. The inflaton and mixed fluctuation spectra 
only appear in the late-time regime, are however sub- 
dominant even there. At early times, t < 20 the Higgs 
spectrum is generated by negative squared masses as in 
tachyonic preheating or quench scenarios. In Ref. [l^ 
it was found that the peak in the momentum spectrum 
pP^^{t,p) can be fitted by a Gaussian; we similarly find, 
at t = 19.5, a spectrum 



with 

A = 0.4482 ±0.0038 (7.42) 

B = 36.1115 ±0.7131 (7.43) 

C = 0.307995 ±0.001168 (7.44) 

At intermediate times the smooth peak broadens and 
decays into spikes, typical of parametric resonance. Para- 
metric resonance also dominates the shape of the spectra 
at late times. As the period and amplitude of oscillation 
change very slowly, the width of the spectrum remains 
constant. 



E. Correlations 

We define the correlation function between the differ- 
ent fluctuations as 



1 



E^i^^ (/r(i,p)/r(t,p)) 



1 



27r2; 



dppsinijir) 



^2^^' (/r(^.p)/r(*'P)) • (7-45) 



We here consider the correlations of the Higgs fluctua- 
tions (i — j — 2), which are displayed in Fig. El We 
observe the correlations to be positive and propagating 
with Ar — 2At, as also found in the large- approxima- 
tion 57]. The propagation with twice the speed of light 
can be related to the fact that the quantum fluctuations 
are correlated by the mean fields whose influence propa- 
gates in opposite space directions. This is corroborated 
by a strong decrease of such correlations when the mean 
fleld amplitude goes to zero jH^ . 



F. Decoherence time 

One of the important questions is the justification of 
using classical instead of quantum dynamics. In the con- 
text of nonequilibrium quantum field theory this has been 
discussed in R.efs . l5 M[60. 61] and applied to the hybrid 
model in Refs. [T^- We use here, adapted to our 
normalization, the definitions of Ref. [ig. The "classi- 
cality" is measured by the imaginary part F{t, p) of a 
correlation function: 



F^j {t, p) = Im 



E 



jr{t,p)ff{t,p) 

2uJa 



(7.46) 



the real part of the bracket being associated with the 
commutator. The criterium for a classical description is 
given by 



pFxJt,p)~^exp [-B{\p\-CY] 



(7.41) 



\F^^\ » 1 



(7.47) 
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FIG. 12: Momentum spectra G,;-, (t, p)p^/(27r^) for the simulation in Fig.0]at the times t — 19.5 (left), t — 51.3 (middle) and 
t = 100.2 (right) 




FIG. 13: Correlation function rC-^xi^i't) for the simulation in 
Fig.H 
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FIG. 14: (Color online) Decoherence time tdGc(p) for which a 
given mode p becomes "classical" (|Jij (tdec, p) | = 1) for the 
simulation in Fig. 2] the red solid line represents the inflaton 
{i ~ j ~ 1) and the green dashed line the Higgs modes (i — 
j = 2); the dotted line corresponds to tdoc oc p^. 



We display in Fig. E|the time for the onset of classicality 
( "decoherence time" ) as a function of momentum. 

As far as the Higgs fluctuations {i = j — 2) are 
concerned the figure can be compared to those of Refs. 
0, 01 . These authors consider the creation of quantum 
fluctuations via the spinodal instability. They simplify 
the time evolution by assuming for the mass of the Higgs 
fluctuations a behavior A^^^ cx (to — t) where Iq marks 
the onset of the spinodal regime, without considering any 



back reactions (see also Refs 0, I^E EE this 
case the mode functions become Airy functions and the 
results are in analytic form. The boundary |F(t,p)| = 1 
between the classical and quantum regimes then behaves 
roughly as t (x . As displayed in Fig. ^^the shape 
of this boundary is quite different in our simulations, the 
classical regimes remains limited within a fixed momen- 
tum band at all times. When one includes back reac- 
tion the behavior of the mass term is linear in time only 
in a very limited time interval; furthermore, due to the 
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inflaton oscillations, the process repeats several times. 
The limited momentum band for which the modes can 
be considered as classical, can be seen as a consequence 
of parametric resonance. Due to the lack of strong dissi- 
pation the oscillations of the classical fields persist at late 
times, and therefore also the resonance band. Whether 
and for which time period this is physical or unphysi- 
cal cannot be determined within the approximation used 
here, even though it is certainly more elaborate as previ- 
ous approaches. 

For the inflaton fluctuations, not considered in Rcf. 9, 
[T^ , the structure of this curve shows that fluctuations at 
very small momenta become classical as early as those of 
the Higgs field, while those for larger momenta develop 
at later times. This is due, presumably, to a stronger 
role of parametric resonance for the evolution of inflaton 
fluctuations. Again the band of momenta for the classical 
regime remains sharply cutoff even at late times. 



VIII. CONCLUSIONS 

In this final section wc would like to comment on some 
of our results and draw some conclusions. 

On the formal level wc have addressed the problem of 
renormalization for a fully coupled system of two quan- 
tum fields in the Hartree approximation, using methods 
similar to those in Refs. Jilj,42j and |5g, but going be- 
yond these applications. We consider this as an essen- 
tial achievement of our work, which will important when 
extending the model by including Goldstone and gauge 
fields and the coupling to gravity where proper renor- 
malization is indispensable [231 ■ Nevertheless we do not 
want to infer that working with judiciously chosen cut- 
offs or less elaborate schemes one may not reach sensible 
physical insights. 

When treated in the Hartree approximation this model 
lacks an efficient mechanism for dissipation. This is con- 
sidered to be a general drawback of the Hartree approx- 
imation. The situation would improve, however, even 
in this approximation, if Goldstone modes were included 
(see, e.g. Fig. 1 in Ref. [l^). Dissipation via particle 
production is also found in the large- iV limit In 
more realistic models dissipation may proceed in addi- 
tion via fermion and gauge fields. The higher order ef- 
fects (sunset graphs or NLO-l/N) become effective only 
when quantum fluctuations have grown to sufficient size. 
This depends of course on the parameters of the model. 
In the earlier stages the Hartree approximation should 
be able to provide reliable information on the evolution 
of decay and resonance processes. The time scales for 
thermalization are expected to be much larger. 

The model displays the expected transition to a broken 
symmetry phase if the initial inffaton amplitude is not 
too large. A possible intermediate restoration of sym- 
metry by quantum ffuctuations would be followed by a 
later transition to the broken symmetry state after fur- 
ther cosmological expansion. Possible consequences like 



the unwanted formation of topological defects have been 
discussed in the hterature |63| . 

We have calculated the boundaries between regions 
where a quantum description is needed and those where 
one may have recurse to classical evolution equations. 
We find marked differences with respect to previous work 
I where the production of quanta is described in a 
simplified way, using a squared mass of the Higgs field 
passing linearly through zero. We find numerically that 
the back reaction limits the classical regime to a low- 
momentum region fixed for all times. Though large exci- 
tations of quantum fluctuations seem to justify the tran- 
sition to a classical description one has to keep in mind, 
that finally one wants to end up with an ensemble de- 
scribed by quantum statistics which is used in the stan- 
dard thermal history of the early universe. The classical 
ensembles suffer from the Rayleigh- Jeans divergence in- 
compatible with the finite amount of initial energy den- 
sity. This automatically forces the fiuctuations back to 
the quantum regime. 

The model discussed here corresponds to the original 
proposal by Linde, using a double well potential. This 
has been used, as a sirnplification or a generic feature in 
some previous studies [3 • For electroweak or GUT scale 
preheating the Higgs sector is based on a symmetry group 
like SU{n) or SO{n) and will in general have more {jih) 
degrees of freedom. Near the spinodal point their masses 
are deg enerate, and this fact has been used in some stud- 
ies jl8| in the way of just using nn identical copies of one 
and the same degree of freedom. However, once the mean 
value of the Higgs field departs from zero, there will be a 
nontrivial mass matrix for the quantum fiuctuations with 
several massless degrees of freedom, the would-be Gold- 
stone bosons. It can be expected that this will modify the 
quantum back reaction in an essential way. While this 
goes beyond the scope of the present investigation, there 
are some studies using classical dynamics with more re- 
alistic Higgs sectors 0, 0, . Our formalism allows 
for a generalization towards more realistic Higgs sectors, 
albeit with the limitation of homogeneous background 
fields. 

To become even more realistic, in particular in view of 
low scale inflation, it would be very desirable to include 
gauge fields as well. Unfortunately any resummation 
or backreaction introduces gauge parameter dependences 
|64i l65l |. which are poorely understood in the framework 
of nonequilibrium quantum field theory. There is some 
recent progress in this direction 66. 67ll68l|. however not 
yet on the level of concrete simulations. 
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APPENDIX A: IDENTITIES FOR FEYNMAN 
INTEGRALS 



An identity that is needed for the renormalization of 
the energy is given by 



P 



(2^) 



D-l 



'"0,a 

327r2 



- In ■ 



''0, a 



(A8) 



Within dimensional regularization {D = A — e) the fol- 
lowing identities (no summation over Greek indices) hold 
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The corresponding Feynman diagrams are depicted in 
Fig. HSl Note that 
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FIG. 15: (a) The Feynman diagram with a topology of a fish 
graph corresponding to Eq. I|A1|I (b) the tadpole type graph 
corresponding to Eq. ljA4^ : in both diagrams the lines de- 
note free propagators with the initial masses mo,ci and mgj^, 
respectively 



APPENDIX B: PERTURB ATIVE EXPANSION 
OF THE MODE FUNCTIONS 

We will present in this section the isolation of the di- 
vergences via a perturbative expansion of the mode func- 
tions (see e.g. Ref. 0i 0) for the case of a coupled 
system of equations. 

Let us split the mode functions /" into a free part 
containing the initial matrix Oij and higher order terms 
represented by the reduced mode functions hf, i.e. 



ftit;p) = e— °*[0,„-)-/i,ni;p)] 



(Bl) 



Here and in the following no summation over Greek in- 
dices is meant if not explicitly stated. 
If we define a potential 

= Ml{t)-MliO) (B2) 

the differential equation (|3.4U|I is equivalent to the fol- 
lowing integral equation 



+ / dt'Klf{t-t';p)VMt')f^it';p) 



fnt;p) 



The retarded kernel of the free equation is given by 



Klf(t-t'-p) = Y,i^Q{t-t')0,pO,p 
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X e 
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Inserting the retarded kernel in the integral equation 
gives 
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where in addition the decomposition in Eq. ljBl|l has been 
used. The dots imply the higher order terms with /i" that 
we do not need for the analysis of the divergences here. 

By partial integration the divergent contributions can 
be isolated in the usual manner 



(B7) 
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This expression and its complex conjugate is all what 
is needed to calculate the divergent contributions in 
Aij{t). We have for the divergent part of the sym- 
metrized Green's function 



2uJa 
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These expressions become divergent if they are integrated 
over (Pp. We will use them as subtraction terms in the 
fluctuation integrals in A[j\t) and E'^^^t) [see Eqs. (14. 6f) 
and gSJ]. 

The divergences in the energy can be found in an anal- 
ogous way by inserting f°'{t,p) and its time derivative in 



the one-loop energy part E'^^^t) [see Eq. (|4.9|) ]. We will 
use an alternative approach in the following. 

The one-loop effective action at time t minus the one 
at t = is given by 

^'Ol (B9) 



□ 



r(i)[M2]^ ^Trln, , 

where it is understood that the numerator and denom- 
inator are 2x2 matrices. This expression can be ex- 
panded locally with respect to V = M^{t) — A^^(O) and 
gradients thereof. The expansion can be obtained by go- 
ing to the momentum representation and by expanding 
with respect to insertions of V(q) and with respect to the 
external momenta q = {qo,q)- As we do not need an in- 
finite wave function renormalization the divergent parts 
are given by the terms of first and second order in V{q). 
We introduce 



(BIO) 

MUt) (Bll) 
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with m^ ^j = A^^^(O). Gq is not the bare propagator 
which would be defined at the vacuum expectation values 
of (j) and X- We diagonalize the initial mass matrix by an 
orthogonal transformation 



or 
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(B12) 
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Then also Gq \j becomes diagonal. The same holds true 
for the inverse matrix. We likewise introduce 



= 0,aG-l{p)0,p , 
where of course 



Vapit))0,p 

(B14) 
(B15) 

is no longer diagonal. The effective action, in the approx- 
imation where all gradient terms are neglected, can now 
be rewritten as 

^Trln{GoG-i} 
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The first terms in the expansion are 
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The three-dimensional reduction is obtained via Eq. (IA4I 
and Eq. (|A1|I . So we find 



p(l) div 



4a;Q 
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The divergent parts of the fluctuation energy are, there- 
fore, 



div 
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a/3 



(B19) 

As a cross check we may obtain the divergent terms in 
the fluctuation integrals which are given by 



A 



(l)dr 



^p(l)div 

Using Eq. (|B15p we have 
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and therefore 
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APPENDIX C: COUNTERTERM STRUCTURE 

In this appendix we will establish a connection between 
the effective counterterms used above and the countert- 
erms in the 2PPI formalism. The latter ones follow from 
a standard counterterm Lagrangian. For details we refer 
the reader to Ref. |2ll5ll53t. 

The general Lagrange density for the model in Eq. (|3.7() 
including a counterterm 5C denotes 



C 



2 ^ 
+5C . 



1 



1 
4! 



The counterterm Lagrangian 5C is given by 

5C = ^^SZijdt'^'d^'^' - ]^5mfj^'¥ 



1 



(C2) 



The gap equation for the renormalized effective mass 
with mass and coupling constant counterterms is 
given by 



+ -{kjki+6kjki){^''^' + . (C3) 



In the 2PPI scheme the following relations for the renor- 
malization constants hold, if one takes a mass indepen- 
dent renormalization scheme 



S'yij-M _ \ij s:/-pq-M 



(C4) 
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Thus all renormalization constants can be derived from 
a vacuum counter term 



5E, 



vac - T^MlMlSC'-'"^' 
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With the help of these identities the gap equation reads 
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with 



Xijkl^n, ~ \jkl^^^ 



2\ijpq5C, 



pq;kl 



'■kl 
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Once the propagator A*^ is renormalized the gap equa- 
tions for the masses A^? become finite. In the Hartree 
approximation, i.e. the one-loop 2PPI approximation, 
the equations for the classical fields 0* are finite as well. 
No wave function renormalization is needed. 

If we speciahze Eq. (jClOp with the help of Eq. (|3.9|l - 
(|3.16l) to our given hybrid model, the renormalization 
constants 5C, simplify to 



x\<t>x ~ ^Cx<t>\x<t> 
\ = ^C<t>x\4>4' 



-2,5C^^ = -2<5exx (Cll) 
-25i^^ (C12) 
. = . (C13) 



(CI) The vacuum counter term in Eq. ljC7p reduces to 



2 ^2 6i^^{Ml^f~25Ux{Ml^f 



(C14) 



Thus the effective counterterm ^A^^ used in Sec. IV CI 
has been mapped to the renormalization constants of a 
standard counterterm Lagrangian 5C [see Eq. (IC2|l ]. In 
principle one can calculate the nonperturbatively fixed 
counterterms 5m^ , 5v^, 5g^ and SX explicitly, insofar 
the renormalized gap equations form a coupled system 
of equations. However, this would not be very enlighten- 
ing. 
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